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Negative-Index Metamaterials: Second-Harmonic Generation, Manley-Rowe 
Relations and Parametric Amplification 
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Second harmonic generation and optical parametric amplification in negative-index metamaterials 
(NIMs) are studied. The opposite directions of the wave vector and the Poynting vector in NIMs 
results in a "backward" phase-matching condition, causing significant changes in the Manley-Rowe 
relations and spatial distributions of the coupled field intensities. It is shown that absorption 
in NIMs can be compensated by backward optical parametric amplification. The possibility of 
distributed-feedback parametric oscillation with no cavity has been demonstrated. The feasibility of 
the generation of entangled pairs of left- and right-handed counter-propagating photons is discussed. 

PACS numbers: 78.67.-n, 42.65.Ky, 42.65.Lm 



I. INTRODUCTION 

Recent demonstration of a negative refractive index 
for metamaterials in the optical range 0, 0] opens new 
avenues for optics and especially nonlinear optics. In par- 
allel with progress for metal-dielectric metamaterials, ex- 
perimental demonstrations of negative refraction in the 
near IR range have been made in a GaAs-based pho- 
tonic crystals Q and in Si-Polyimide photonic crystals 
0. Negative refractive- index metamaterials (NIMs) are 
also referred to as left-handed materials (LHMs). The 
sufficient (but not necessary) condition for a negative re- 
fractive index is simultaneously negative dielectric per- 
mittivity e(u>) and negative magnetic permeability [i(uj) 
0. Negative magnetic permeability in the optical range 
has been demonstrated in 

H0||. NIMs exhibit highly 
unusual optical properties and promise a great variety 
of unprecedented applications. Optical magnetization, 
which is normally ignored in linear and nonlinear-optics 
of the ordinary, positive-index materials (PIMs) (i.e., 
right-handed materials, RHMs) plays a crucial role in 
NIMs. 

The main emphasis in the studies of NIMs has been 
placed so far on linear optical effects. Recently it has 
been shown that NIMs including structural elements with 
non-symmetric current-voltage characteristics can pos- 
sess a nonlinear magnetic response at optical frequencies 
H and thus combine unprecedented linear and nonlinear 
electromagnetic properties. Important properties of sec- 
ond harmonic generation (SHG) in NIMs in the constant- 
pump approximation were discussed in |lOl | for semi- 
infinite materials and in [llj for a slab of a finite thick- 
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ness. The propagation of microwave radiation in non- 
linear transmission lines, which are the one-dimensional 
analog of NIMs, was investigated in The possi- 

bility of exact phase-matching for waves with counter- 
propagating energy-flows has been shown in |l3l | for the 
case when the fundamental wave falls in the negative- 
index frequency domain and the SH wave lies in the 
positive-index domain. The possibility of the existence of 
multistable nonlinear effects in SHG was also predicted 
in 0. 

As seen from our consideration below, the phase- 
matching of normal and backward waves is inherent for 
nonlinear optics of NIMs. We note here that the impor- 
tant advantages of interaction schemes involving counter- 
directed Poynting vectors in the process of optical para- 
metric amplification in ordinary RHMs were discussed in 
early papers ^4[. However, in RHMs such schemes im- 
pose severe limitations on the frequencies of the coupled 
waves because of the requirement that one of the waves 
has to be in the far-infrared range. 

Absorption is one of the biggest problems to be ad- 
dressed for the practical applications of NIMs. In [Tol[lll| . 
a transfer of the near-field image into SH frequency do- 
main, where absorption is typically much less, was pro- 
posed as a means to overcome dissipative losses and thus 
enable the superlens. 

In this paper, we demonstrate unusual characteristics 
in the spatial distribution of the energy exchange between 
the fundamental and second- harmonic waves. Both semi- 
infinite and finite-length NIMs are considered and com- 
pared with each other and with ordinary PIMs. Our anal- 
ysis is based on the solution to equations for the coupled 
waves propagating in lossless NIMs beyond the constant- 
pump approximation. The Manley-Rowe relations for 
NIMs are analyzed and they are shown to be strikingly 
different from those in PIMs. We also propose a new 
means of compensating losses in NIMs by employing opti- 
cal parametric amplification (OPA). This can be realized 
by using control electromagnetic waves (with frequencies 
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outside the negative- index domain), which provide the 
loss-balancing OPA inside the negative-index frequency 
domain. We also predict laser oscillations without a cav- 
ity for frequencies in the negative-index domain and the 
possibility of the generation of entangled pairs of counter- 
propagating right- and left-handed photons. 

The paper is organized as follows. Section HT1 discusses 
the unusual spatial distribution of the field intensities 
for SHG in finite and semi- infinite slabs of NIMs. The 
Manlcy-Rowc relations are derived and discussed here. 
The feasibility of compensating losses in NIMs by using 
the OPA is considered in Section ITTT1 In this Section we 
also study cavity-less oscillations based on distributed 
feedback. Finally, a summary of the obtained results 
concludes the paper. 



II. SECOND HARMONIC GENERATION IN 
NIMS 

A. Wave vectors and Poynting vectors in NIMs 

We consider a loss-free material, which is left-handed 
at the fundamental frequency to\ (t\ < 0, /.ti < 0), 
whereas it is right-handed at the SH frequency cu 2 = 2ui 
(t 2 > 0, fi2 > 0). The relations between the vectors of 
the electrical, E, and magnetic, H, field components and 
the wave- vector k for a traveling electromagnetic wave, 

E(r,i) = E (r)expH(orf- k- r)] +c.c, (1) 
H(r,i) = H (r)exp [-i (art-k- r)] +c.c, (2) 



are given by the following formulas 



k x E 

V~eE{r,t) 



(w/c)/xH, k x H 
-V^ff(r, t), 



-(w/c)eE, (3) 
(4) 



which follow from Maxwell's equations. These expres- 
sions show that the vector triplet E, H and k forms a 
right-handed system for the SH wave and a left-handed 
system for the fundamental beam. Simultaneously neg- 
ative ti < and fa < result in a negative refractive 
index n = —^[Jit. As seen from Eqs. Q and J5J), the 
phase velocity ~v p h is co-directed with k and is given by 
v ph = (k/AOjw/fc) = (k/fc)(c/M), where k 2 = ii 2 (lu/c) 2 . 
In contrast, the direction of the energy flow (Poynting 
vector) S with respect to k depends on the signs of e and 
fj,; 



S(r,t) 



-[Ex H] = 

Air Anujt 



c 2 k 

AlTLUt 



[H x k x H] 



c 2 k 



-E z 



(5) 



As mentioned, we assume here that all indices of e, \x and 
n are real numbers. Thus, the energy flow Si at cJi is 
directed opposite to k 1; whereas S2 is co-directed with 
k 2 . 



B. SHG: Basic equations and the Manley-Rowe 
relations 



We assume that an incident flow of fundamental ra- 
diation Si at u>i propagates along the z-axis, which is 
normal to the surface of a metamaterial. According to 
(JSJ, the phase of the wave at u>\ travels in the reverse 
direction inside the NIM (the upper part of Fig|TJ . Be- 
cause of the phase-matching requirement, the generated 
SH radiation also travels backward with energy flow in 
the same backward direction. This is in contrast with the 
standard coupling geometry in a PIM (the lower part of 
FigEJ. Following the method of 0], we assume that a 
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FIG. 1: SHG geometry and the difference between SHG in 
LHM and RHM slabs. ' 

nonlinear response is primarily associated with the mag- 
netic component of the waves. Then the equations for 
the coupled fields inside a NIM in the approximation of 
slow- varying amplitudes acquire the form: 



dA 2 
dz 

dAi 
dz 

P) 



*?4 4 ^f / ) / ^ cx p(- Afcz )' 



k\c 2 



torgjjA^Al exp(Afcz) 



(6) 



Here, X e /f is the effective nonlinear susceptibility, Ak = 
k 2 — 2ki is the phase mismatch, and A 2 and A\ are the 
slowly varying amplitudes of the waves with the phases 
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traveling against the z-axis: 

Hj(z, t) = Aj exp[—i(kjZ + u>jt)] + c.c. 



(8) 



where, u> 2 — 2wi and k\ t2 > are the moduli of the wave- 
vectors directed against the z-axis. We note that accord- 
ing to Eq. (0} the corresponding equations for the electric 
components can be written in a similar form, with tj sub- 
stituted by Hj and vice versa. The factors fij were usually 
assumed to be equal to one in similar equations for PIMs. 
However, this assumption does not hold for the case of 
NIMs, and this fact dramatically changes many conven- 
tional electromagnetic relations. The Manley-Rowe rela- 
tions for the field intensities and for the energy flows 
follow from Eqs. © - 0: 



hdlA^ 2 k 2 d\A 2 {< 



£\ dz 2&2 dz 



= 0, 



d|^| 2 d\S 2 f 



dz 



dz 



= 0. (9) 



The latter equation accounts for the difference in the 
signs of ei and e 2 , which brings radical changes to the 
spatial dependence of the field intensities discussed be- 
low. 

We focus on the basic features of the process and ig- 
nore the dissipation of both waves inside the nonlinear 
medium; in addition, we assume that the phase match- 
ing condition k 2 = 2ki is fulfilled. The spatially-invariant 
form of the Manley-Rowe relations follows from equation 
(23): 



|A 1 | 2 Ai + l^ 2 | 2 A2 = c, 



(10) 



where C is an integration constant. With e\ = —e 2 , 
which is required for the phase matching, equation H10|) 
predicts that the difference between the squared ampli- 
tudes remains constant through the sample 



c. 



(11) 



as schematically depicted in the upper part of Fig. ^ 
This is in striking difference with the requirement that 
the sum of the squared amplitudes is constant in the 
analogous case in a PIM, as schematically shown in the 
lower part of Fig. ^ We introduce now the real phases 
and amplitudes as Ai 2 = h\ i2 cxp(i(/>i.2). Then the equa- 
tions for the phases, which follow from Eqs. © and J7J, 
show that if any of the fields becomes zero at any point, 
the integral I|1U|) corresponds to the solution with the 
constant phase difference 2<^>i — <p 2 = tt/2 over the entire 
sample. 

The equations for the slowly-varying amplitudes cor- 
responding to the ordinary coupling scheme in a PIM, 
shown in the lower part of Fig. ^ arc readily obtained 
from Eqs. © - JSJ by changing the signs of k\ and 
k 2 . This does not change the integral (|10[) : more im- 
portantly, the relation between ei and e 2 required by the 
phase matching now changes to ei = e 2 , where both con- 
stants arc positive. The phase difference remains the 
same. Because of the boundary conditions h\(0) = h±o 
and h 2 (0) = Ii 2 q = 0, the integration constant becomes 



C = h 2 . Thus, the equations for the real amplitudes in 
the case of a PIM acquire the form: 



hi(z) = ^fi 2 w - h 2 {z) 2 , 
dh 2 /dz = ft[/i 2 - h 2 (z) 2 ] 



(12) 
(13) 



with the known solution 



h 2 (z) = h w tanh(z/z ), (14) 
hi(z) = h 10 / cosh(z / z ) , z = [k/ii ] _1 . (15) 

Here, k = (e 2 (jj 2 / k 2 c 2 )AirXej t ■ The solution has the same 
form for an arbitrary slab thickness, as shown schemati- 
cally in the lower part of Fig. ^ 



C. SHG in a NIM slab 

Now consider phase-matched SHG in a lossless NIM 
slab of a finite length L. Equations © and 111|) take the 
form: 



h^z) 2 = C + h 2 {z) 2 , 
dh 2 /dz = -k[C + h 2 (zf] 



(16) 
(17) 



Taking into account the different boundary conditions in 
a NIM as compared to a PIM, hi(0) = /iio an d h 2 (L) = 0, 
the solution to these equations is as follows 



(18) 
(19) 



h 2 = VCtan[vC«;(i- z)], 
hi = VC/cos[VCk(L- z)], 



where the integration parameter C depends on the slab 
thickness L and on the amplitude of the incident funda- 
mental radiation as 



(20) 



CkL = cos^iVC/hw). 



Thus, the spatially invariant field intensity difference be- 
tween the fundamental and SH waves in NIMs depends 
on the slab thickness, which is in strict contrast with the 
case in PIMs. As seen from equation l|39|) . the integra- 
tion parameter C = hi(z) 2 — h 2 (z) 2 now represents the 
deviation of the conversion efficiency r\ = h 2Q /h 2 from 
unity: (C/h\ Q ) = 1— 77. Figure|2]shows the dependence of 
this parameter on the conversion length zq = (nhio)^ 1 . 
The figure shows that for the conversion length of 2.5, 
the NIM slab, which acts as nonlinear mirror, provides 
about 80% conversion of the fundamental beam into a 
reflected SH wave. Figure [21 depicts the field distribution 
along the slab. One can see from the figure that with an 
increase in slab length (or intensity of the fundamental 
wave) , the gap between the two plots decreases while the 
conversion efficiency increases (comparing the main plot 
and the inset). 
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FIG. 2: The normalized integration constant C/h\ and the 
energy conversion efficiency 77 vs the normalized length of a 
NIM slab. 




FIG. 3: The squared amplitudes for the fundamental wave 
(the dashed line) and SHG (the solid line) in a lossless NIM 
slab of a finite length. Inset: the slab has a length equal 
to one conversion length. Main plot: the slab has a length 
equal to five conversion lengths. The dash-dot lines show the 
energy-conversion for a semi-infinite NIM. 



D. SHG in a semi-infinite NIM 

Now we consider the case of a semi-infinite NIM at 
z > 0. Since both waves disappear at z — > 00 due to 
the entire conversion of the fundamental beam into SH, 
C = 0. Then equations l|39|) and l|17|) for the amplitudes 
take the simple form 

h 2 (z) = h 1 (z), (21) 
dh 2 /dz = -nh\. (22) 

Equation (|21|1 indicates 100% conversion of the incident 
fundamental wave into the reflected second harmonic at 
z = in a lossless semi-infinite medium provided that 
the phase matching condition Afc ~ is fulfilled. The 



integration of with the boundary condition hi (0) = 
h\o yields 

W = 7 7 ^— ,z = ( K h w )- 1 . (23) 
{z/zo) + 1 

Equation (|23|l describes a concurrent decrease of both 
waves of equal amplitudes along the z-axis; this is shown 
by the dash-dot plots in Fig. For z ^> z , the depen- 
dence is inversely proportional to z. These spatial depen- 
dencies, shown in Fig. [3J are in strict contrast with those 
for the conventional process of SHG in a PIM, which are 
known from various textbooks (compare, for example, 
with the lower part of Fig^l . 

III. OPTICAL PARAMETRIC AMPLIFICATION 
AND DIFFERENCE-FREQUENCY 
GENERATION IN A NIM SLAB WITH 
ABSORPTION 

A. OPA: basic equations and Manley-Rowe 
relations 

As mentioned in Subsection III Al S is countcrdirected 
with respect to k in NIMs, because e < and [i < 0. 
We assume that a left-handed wave at uji travels with its 
wave-vector directed along the z-axis. Then its energy 
flow Si is directed against the z-axis. We also assume 
that the sample is illuminated by a higher-frequency 
electromagnetic wave traveling along the axis z. The 
frequency of this radiation lu 3 falls in a positive index 
range. The two coupled waves with co-directed wave- 
vectors k3 and ki generate a difference- frequency idler at 
LU2 = UJ3 — lui, which has a positive refractive index. The 
idle wave contributes back into the wave at uj\ through 
three-wave coupling and thus enables optical paramet- 
ric amplification (OPA) at uj\ by converting the energy 
of the pump field at W3. Thus, the nonlinear-optical pro- 
cess under consideration involves three- wave mixing with 
wave- vectors co-directed along z. Note that the energy 
flow of the signal wave, Si, is directed against z, i.e., it 
is directed opposite to the energy flows of the two other 
waves, S2 and S3 (Fig. 0] the left part). Such a coupling 
scheme is in contrast with the ordinary phase-matching 
scheme for OPA, which is schematically shown in the 
right part of Fig. 21 As above, we consider the magnetic 
type of the quadratic nonlincarity. For the magnetic field 

Hj(z, t) = hj cxp[i(kjZ — ujjt)] + c.c, (24) 

the nonlinear magnetization at the signal and idler fre- 
quencies is given by the equations 

M? L = 2 X % f h 3 h* 2 exp{i[(k3~k 2 )-u 1 t}} 1 (25) 

M" L = 2 X ^ f h 3 ht exp{i[(fc 3 -fci)-w 2 t]}. (26) 

Here, j = 1, 2, 3; u> 2 = u>3 — toi; and kj — \rij\wj/c > 0. 
Then the equations for the slowly-varying amplitudes of 
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FIG. 4: The difference between OPA processes in LHM and 
PIM slabs. 



the signal and the idler acquire the form 

dhi . „ r . . a\ 
= Kj\hj,h 2 exp iAkz H 

GtZ 2 

* 2 i j # r . A , , a 2 



where <7j = 8irxfff e j LL >'j/kjC 2 , Ak 



A:, 



(27) 
(28) 

k 2 — ki, and 



ctj are the absorption indices. The amplitude of the 
pump /i 3 is assumed constant. We note the following 
three fundamental differences in equation l|27[) as com- 
pared with the ordinary difference-frequency generation 
(DFG) through the three-wave mixing of co-propagating 
waves in a PIM. First, the sign of o\ is opposite to that 
of <r 2 because t\ < 0. Second, the opposite sign appears 
with ai because the energy flow Si is directed against 
the z-axis. Third, the boundary conditions for h\ are 
defined at the opposite side of the sample as compared 
to h 2 and /13 because their energy-flows Si and S2 are 
counter-directed. 

At ct\ = a-2 = 0, one finds with the aid of Eqs. (I27II . 
and JSJ: 



d_ 

dz 
d 
dz 



Slz S22 

hwi hijji 

ei wi 



0, 

\h-2 



£2 ^2 



= 0. 



(29) 
(30) 



These equations represent the Manley-Rowe relations 
|l5| . which describe the creation of pairs of entangled 
counter-propagating photons f\w\ and hwi- The equa- 
tions account for the opposite sign of the corresponding 
derivatives with respect to z. Equation (|3U|) predicts that 
the sum of the terms proportional to the squared ampli- 
tudes of signal and idler remains constant through the 
sample, which is in contrast with the requirement that 
the difference of such terms is constant in the analogous 
case in a PIM. We note that according to Eqs. © and 
|JHJ the corresponding equations for the electric compo- 
nents in the case of the quadratic electric nonlincarity can 
be written in a similar form with ej substituted by fij. 



As seen from the equations below, this does not change 
either the results obtained or the main conclusions pre- 
sented here; the same is true for the case of SHG. As 
mentioned in Section^ the factors pLj were usually as- 
sumed equal to unity in equations for PIMs, which is not 
the case for NIMs. 



B. OPA and DFG in NIMs 



We introduce the normalized amplitudes a,j = 
\Jtjl fXjhj/^/uJJ; their squared values are proportional to 
the number of photons at the corresponding frequencies. 
The corresponding equations for such amplitudes acquire 
the form 



dai . „ m 

—j— = — iga 2 exp[iAfczJ + — a±, 

da 2 . „ r . A , a 2 
= iga x exp[iAkz\ — — a 2 , 



(31) 
(32) 



where g = (^/ZZJIZJi"/ \J 'eie 2 / '(^ I 'c)x^ 2) ^3- Account- 
ing for the boundary conditions a\{z = L) = an,, and 
a 2 (z = 0) = a2o (where L is the slab thickness), the 
solutions to equations (|31|l and l|32|) are as follows 

Ak 

ai(z) = A 1 exp[(/3i + i— )z] + 
Ak 

+ A 2 exp[(/3 2 + i—)z], (33) 
Afc 

a* 2 (z) = KiAiexp[(J3i—i—)z] + 

Ak 

+ K 2 A 2 exp[(/3 2 -i— )z]. (34) 

Here. 

[3 h2 = («! - a 2 )/4 ± iR, Ki,a = [±R + is]/g, (35) 

R = vV - s 2 , s = (ai + a 2 )/4 - iAk/2, (36) 

Afc 

Ai = {a 1L n 2 - a* 20 exp[(/3 2 + i—)L]}/D, (37) 

Afc 

A 2 = -{a 1L K! - a* 20 exp[(/3i + i—)L]}/D, (38) 

Afc Afc 
D = K 2 exp[(A + i—)L] - m exp[(/3 2 + i— )L}. (39) 

At a2o = 0, the amplification factor for the left-handed 
wave is given by rj a (ui\) = |aio/aii| 2 , where 



aw 
aiL 



exp 




f ai — ol 2 Afc 

V % — 

x 4 2 


y. 




cos RL + 


ol\ + oi2 . Ak 
4R l 2R 


sin RL 



(40) 



Alternatively, at a* 1L = 0, the conversion factor for the 
difference-frequency generation of the left-handed wave is 
found as r) g (u>\) = |aio / Q-Iq | 2 , where 



aio 

a 20 



-(g/R) sin RL 



cos RL 



ai + a 2 .Afc 
4i? % 2R 



(41) 



sin RL 
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Equation i|40|) shows that the amplification of the left- 
handed wave can be turned into a cavity-less oscillation 
when the denominator tends to zero. The conversion 
factor for DFG, r/ gi experiences a similar increase. In the 
case of Afc = and small optical losses (a± + a 2 )L <C 7r, 
equations (pT3|> and are reduced to 



at(z) 
a 2 (z) 



'IL 



cos(gL) 
ia* 



"XL 



cos(gz) 
s'm(gz) 



^20 

cos(gL) 

O20 



cos(gL) cos(gL) 
The output amplitudes are then given by 



sin[g(z 
cos[g(z 



L 1L 



no 



0-2L 



cos(gL) 
ia* L tan(gL) 



ia 2 otan(gL), 
a-20 



cos(gL) 



mm 

L){43) 

(44) 
(45) 



Thus, the oscillation threshold value for the control field 
intensity in this case is given by g t — ir/2L. It increases 
with absorption and phase mismatch. 
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FIG. 6: Resonant changes in the distribution of the normal- 
ized intensity of the left-handed wave inside the slab of NIM, 
rja(uji) (the solid line) and n g (uJi) (the dashed line), caused 
by the adjustment of the normalized intensity for the control 
field at lu 3 , gL. a x L = 1, a 2 L = 1/2, Ak = 0. 
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FIG. 5: The amplification factor r/ a (uji) (the solid line) and 
the efficiency of difference-frequency generation r/ g (uii) (the 
dashed line) for the backward wave at z~0. ol\L = 1, otiL = 
1/2. The upper plot: Ak = 0. The lower plot Ak = n. 

The dependence of the output intensity for the left- 
handed wave propagating in an absorptive NIM slab in 
the presence of the control field at W3 and at 020 = is 
shown with the solid line in Fig. [5] for two representa- 
tive cases of exact and partial phase- matching. The dash 
plots show the output in the case of DFG (at oil = 0, 
a 2,a =/= 0). Amplification in the upper part of Fig. [S] 



reaches many orders in the first maximum and increases 
in the next maximums. It is seen that the amplification 
can entirely compensate for absorption and even turn into 
oscillations when the intensity of the control field reaches 
values given by a periodic set of increasing numbers. The 
larger the corresponding value, the greater is the amplifi- 
cation and the DFG output; the latter depends on the ab- 
sorption for both waves and on the phase mismatch Afc. 
The conversion factor is larger in its maximums than the 
amplification factor because DFG is a one-step process, 
whereas OPA is a two-step process as discussed in Sub- 
section 1^^] As seen from Fig. the output shows a 
resonance dependence on the intensity of the control field 
at W3. Figure El depicts corresponding changes in the dis- 
tribution of the negative-index field inside the slab in the 
vicinity of the first maximum at Afc = 0. 



IV. CONCLUSION 

We have studied the unusual properties of second- 
harmonic generation (SHG) in mctamatcrials that have 
a negative refractive index for the fundamental wave and 
a positive index for its second harmonic (SH). The pos- 
sibility of a left-handed nonlinear-optical mirror, which 
converts the incoming radiation into a reflected beam 
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at the doubled frequency with efficiency that can ap- 
proach 100% for lossless and phase- matched medium is 
considered. The most striking differences in the nonlinear 
propagation and the spatial dependence of the energy- 
conversion process for SHG in NIMs, as compared to 
PIMs, can be summarized as follows. In NIMs, the in- 
tensities of the fundamental and SH waves both decrease 
along the medium. Such unusual dependence and the 
apparent contradiction with the ordinary Manlcy-Rowc 
relations are explained by the fact that the energy flows 
for the fundamental and SH waves are counter-directed, 
whereas their wave-vectors are co-directed. Another in- 
teresting characteristic of SHG in NIMs is that the energy 
conversion at any point within a NIM slab depends on 
the total thickness of the slab. This is because SHG in a 
NIM is determined by the boundary condition for SH at 
the rear interface rather than the front interface of the 
slab. 

We have shown the feasibility of compensating losses in 
NIMs by optical parametric amplification (OPA). In this 
process, the wave-vectors of all three coupled waves are 
co-directed, whereas the energy flow for the signal wave is 
counter-directed with respect to those for the pump and 
the idler waves. This process is characterized by prop- 
erties that are in strict contrast with those known for 
conventional nonlinear-optical crystals. Such extraordi- 
nary features allow one to realize optical parametric os- 
cillations (OPOs) without a cavity at frequencies where 
the refractive index is negative. We also showed that the 



OPA and OPO in NIMs enable the generation of pairs 
of entangled counter-propagating right- and left-handed 
photons inside the NIM slabs. 

The backward energy flow for one of the coupled waves 
(whereas the wave- vectors of all the coupled waves are 
co-directed) is inherent for NIMs and it makes this pro- 
cess different from three-wave mixing in PIMs. This is 
also different from various processes in RHMs based on 
distributed gratings and feedback. The important ad- 
vantage of the backward OPA and OPO in NIMs inves- 
tigated here is the distributed feedback, which enables 
oscillations without a cavity. In NIMs, each spatial point 
serves as a source for the generated wave in the reflected 
direction, whereas the phase velocities of all the three 
coupled waves are co-directed. As mentioned, it is very 
hard to realize such a scheme in PIMs, while the OPA in 
NIMs proposed herein is free from the PIM limitations. 
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